Phonon-Induced Spin Relaxation of Conduction Electrons in Aluminum 



Jaroslav Fabian and S. Das Sarma 
Department of Physics, University of Maryland at College Park, College Park, Maryland 20742-4111 

Spin-flip Eliashberg function a%F and temperature-dependent spin relaxation time Ti(T) are calcu- 
lated for aluminum using realistic pseudopotentials. The spin-flip electron-phonon coupling constant 
As is found to be 2.5 x 10~^. The calculations agree with experiments validating the Elliott- Yafet 
theory and the spin-hot-spot picture of spin relaxation for polyvalent metals. 
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PACS numbers: 71.70.Ej, 75.40. Gb, 76.30. Pk 

Spin dynamics of itinerant electrons in metals and 
semiconductors is attracting increasing attention. Part 
of the reason for this interest is fundamental, arisin, 
from improved spin injection and detection techniques 
which now allow precise measurements of spin transport, 
relaxation, and coherence properties. But much of the 
recent interest is also motivated by the exciting potential 
of using electron spin as a building block in nanoelec- 
tronics (dubbed "spintronics" ) where spin dynamics and 
transport is projected to be utilized in proposed novel 
device applications. The most ambitious such possibility 
is using electron spin as a qubit in a quantum computer 
architecture, but more modest proposals involving the 
use of spin injection and transport in new quantum tran- 
sistor devices ("spin transistors") have also been made 

0. 

Electron spin already plays a fundamental, albeit pas- 
sive, role in giant magnetoresistance-based memory de- 
vices. The current push for a better understanding of 
spin dynamics in electronic materials is, however, based 
on the hope that the electron spin could be used as an 
active element, where manipulation of spin in a con- 
trolled manner will lead to novel device applications 
which are not feasible in conventional microelectronics. 
This hope arises from two underlying concepts: the inher- 
ently quantum mechanical nature of spin (enabling the 
possibility of truly quantum devices which could not be 
envisioned within standard micro- or nanoelectronics), 
and, even more importantly, the inherently long relax- 
ation or coherence time of spin eigenstates in metals and 
semiconductors (indeed, in a typical nonmagnetic metal 
at room temperature electron spins survive for hundreds 
of picoseconds; by comparison, momentum states live 
no more than femtoseconds). This Letter provides the 
first realistic quantitative calculation of the temperature 
dependent spin relaxation time (the so called Ti relax- 
ation time) in an electronic material, namely, metallic 
aluminum. The calculation, for reasons to be explained 
below, is surprisingly subtle and extremely computation- 
ally demanding; it has therefore never been attempted 
before although the basic theory for the phenomenon 
goes back more than thirty- five years H,^. 

The mechanism behind spin relaxation in metals is 
believed to be the spin-flip scattering of electrons off 
phonons and impurities, as suggested by Elliott j^] and 



Yafet The periodic, ion-induced spin-orbit interac- 
tion causes electronic Bloch states to have both spin 
up and spin down amplitudes. The states can still 
be polarized by a magnetic field (so we can call them 
up and down) but because of the spin mixing, even a 
spin-independent interaction with phonons or impurities 
(which are assumed to be nonmagnetic) leads to a tran- 
sition from, say, up to down, degrading any unbalanced 
spin population. (Note that the spin-orbit interaction by 
itself does not produce spin relaxation-what is needed is 
spin-orbit coupling to mix the up and down spins, and 
a momentum conservation-breaking mechanism such as 
impurities or phonons.) Although these arguments seem 
to be consistent with experimental findings, there has 
been to date no calculation of Ti for a metal based on 
the Elliott- Yafet theory. 

In this Letter we calculate the phonon contribution to 
Ti for aluminum providing the first quantitative justifica- 
tion of the theory. (Impurities in real samples contribute 
only a temperature independent background which can 
be subtracted from the measurement.) At tempera- 
tures T above 100 K, where experimental data are not 
available, our calculation is a prediction which should 
be useful for designing room-temperature spintronic de- 
vices that use aluminum. We also calculate the spin-flip 
Eliashberg function Q;|.F(il) which measures the ability 
of phonons with frequency Q to change electron momenta 
and spins. This function, which is an analogue of the or- 
dinary (spin-conserving) Eliashberg function a^F{Q) [Q, 
is important in spin-resolved point-contact spectroscopy 
where phonon-induced spin flips could be directly ob- 
served. (A recent effort to detect phonon-induced 
spin flips in aluminum failed because of the overwhelming 
spin-flip boundary scattering in the sample.) 

Aluminum belongs to the group of metals whose 
spin relaxation is strongly influenced by band-structure 
anomalies . Monod and Beuncu Q observed that while 
simple estimates based on the Elliott- Yafet theory work 
well for monovalent alkali and noble metals, they severely 
underestimate l/Ti for polyvalent Al, Mg, Be, and Pd 
(the only polyvalent metals measured so far). Silsbee 
and Beuneu |^ pointed out that in aluminum accidental 
degeneracies can significantly enhance 1 /Ti . We recently 
Q developed a general theory including band structure 
anomalies like accidental degeneracy, crossing Brillouin 
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zone boundaries or special symmetry points, and rigor- 
ously showed that they all enhance 1/Ti. This explains 
the Monod-Beuneu finding because the anomalies (which 
we named "spin hot spots" |^) are ubiquitous in poly- 
valent metals. The present calculation is consistent with 
the spin-hot-spot picture. 

The formula for the spin relaxation rate, first derived 
by Yafet Q, can be written in the more conventional 
electron-phonon terminology as 

l/T,{T) = %^T dnalF{n)^^^ (1) 

where N{n) = [exp(?irj/fcsT) - 1]-^ and is 
the spin-fiip Eliashberg function. Before writing the ex- 
pression for OigF we introduce the following notation. 
Electron states 4* (normalized to a primitive cell) in the 
periodic potential V containing the spin-orbit interac- 
tion are labeled by lattice momentum k, band index n, 
and spin polarization | or |. If has inversion sym- 
metry (as in aluminum), states VPknT s-nd ^kni are de- 
generate The spin polarization then means that 
these two states are chosen to satisfy (vl/knt; '^'z^knt) = 
— (^kni, o'z^'knj,) > with the off-diagonal terms van- 
ishing p] . Lattice vibrations are represented by phonons 
with momentum q and polarization index u. Phonon fre- 
quency is and polarization vector Uq^ (we consider 
a Bravais lattice). If q = k — k' and 

5k„T,k'n'i = |Uq. • (*k„T, Vy*k'„'i)l' : (2) 

the spin-flip Eliashberg function is 

4^(f^) - ^ E«5k„T,k'„'i'5(c^q,. - f^))k„)k'„'. (3) 

V 

Here gs is the number of states per spin and atom at the 
Fermi level, M is the ion mass, and (...}kn denotes the 
Fermi surface averaging jl^ . 

We calculate o?gF and T\ for aluminum by the pseu- 
dopotential method The spin- independent part of 
the electron-ion pseudopotential is represented by the 
Masovic-Zekovic pT| semi-empirical form factor which 
reproduces well the observed band gaps at the symme- 
try points of the Brillouin zone. This is a crucial fea- 
ture because the presence of spin hot spots makes T\ 
sensitive to the band structure at the Fermi surface 
The spin-orbit part of the pseudopotential comes from a 
fit of the first-principles Bachelet-Hamann-Schliiter pseu- 
dopotential ||l2| to ckL • ST^i, where L (S) is the orbital 
(spin) momentum operator and Vi is the operator pro- 
jecting on the orbital momentum state /. The parameter 
a — 2A X 10~^ a. u. (1 a. u. =2 Ry) inside the ion 
core of twice the Bohr radius, rc = 2rB- Outside the 
core a vanishes. The cutoff for the plane- wave energy is 
1 a. u. from the Fermi level fl^ . For phonons we use 
the highly successful force-constant model of Cowley 
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FIG. 1. Calculated spin relaxation time Ti of conduction 
electrons in aluminum as a function of temperature T (solid 
line). Symbols are experimentally determined phonon 
contribution to Ti from measurements by Johnson and Silsbee 
[ p^ (triangles) and Lubzens and Schultz (squares). The 
dashed line is an estimate of Ti from Eq. 07 The inset shows 
Ti over a wider temperature range with thin lines indicat- 
ing the predicted low-T (Ti ~ T"^) and high-T (Ti ~ T"^) 
behavior. 

which gives an excellent fit to the experimental spectrum. 
Finally, the sums over the Brillouin zone are done by the 
tetrahedron method with a specially designed grid 
of more than 4000 points around the Fermi surface in 
an irreducible wedge of the Brillouin zone to accurately 
obtain contributions from the spin hot spots. 

Figure ^ shows the calculated spin relaxation time 
Ti as a function of temperature. The agreement with 
experiment is evident. At high temperatures where 
there are no experimental data, our calculation predicts 
ri[ns] « 2iT~^[K]. This behavior is expected for a 
phonon-induced relaxation above the Debye temperature 
which for aluminum is about 400 K. As Fig. |l| shows the 
Ti ~ behavior starts already at 200 K. At very low 
temperatures the theory predicts the asymptotic tem- 
perature dependence Ti T^^ (the Yafet law |^) purely 
on dimensional grounds. Our calculation gives rather a 
good fit to Ti ~ T-^-^^ between 2 and 10 K. At lower 
temperatures our results cease to be reliable because of 
the finite size (limited by the computing resources) of the 
tetrahedron blocks in the summations over the Brillouin 
zone. We anticipate that the asymptotic Yafet law would 
be reached at lower temperatures (much lower than 2 K) 
since we have verified numerically its origin, namely that 
.9k'nT,k'ni ~ (k - k')"^ at k ^ k' 1^ (a quadratic de- 
pendence would be expected for spin-conserving matrix 
elements). In Fig. ^we also plot an estimate of Ti based 
on the simple formula |^ 

Ti«t/4(62}, (4) 
where {b"^) is the Fermi surface average of the spin-mixing 
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parameter, calculated in Q to be 2.5 x 10~^, and r is 
the momentum relaxation time obtained from the Drude 
formula for the resistivity (resistivity data taken from 
Ref. ||l8|) with an electron thermal mass of 1.5 of 
the free electron mass. This estimate of Ti reproduces 
well the calculated functional temperature dependence 
making Eq. ^ useful as a starting point for order-of- 
magnitude estimates. 

The calculated spin-flip Eliashberg function oigF for 
aluminum is shown in Fig. ^ along with the phonon 
density of states F and the spin-conserving Eliashberg 
function a^F. The last agrees very well with previous 
calculations |^,^. Transverse phonon modes which 
dominate the low-frequency spectrum are less effective 
in scattering electrons, with or without spin flip, than 
high-frequency longitudinal phonon modes. The behav- 
ior of a^F at small ft that gives the Yafet law is predicted 
to be OgF ^ ri^. We are not able to reproduce this re- 
sult, again because of the finite size of the tetrahedron 
blocks. This is a well known problem that the asymptotic 
low- frequency behavior is hard to reproduce (2^j21 ] . 

From the Eliashberg function we can calculate the ef- 
fective electron-phonon coupling constant 



A(5) = 2 



dQ 2 



F{n). 



(5) 



We obtain A w 0.4 and As w 2.5 x 10"^ The spin- 
conserving A falls well into the interval of the "rec- 
ommended" values 0.38 — 0.48 obtained by differ- 
ent methods p0| , pl| , p2|j2^ . At high temperatures the 
phonon-induced relaxation is determined by A(5), since 
in this regime h/r « 27rAfcsr and h/Ti « AnXskBT. 
The momentum to spin relaxation time ratio t/Ti is 
2A5/A w 1.24 x 10""*. From the above ratio of r/Ti 
we obtain the "effective" (6^) « 3.1 x 10~^ in Eq. not 
that different from its calculated value of 2.0 x 10^ 
Thus, our theory is internally consistent. 

We conclude with a remark on the accuracy of our 
calculation of Ag. The numerical error is accumulated 
mostly during the summations over the Brillouin zone. 
This error was previously estimated Q to be about 
10%. Another source of uncertainty, which is much 
more important here than in the spin- conserving calcu- 
lations, comes from the choice of the pseudopotentials. 
While the spin-orbit pseudopotential sets the overall scale 
(1/Ti ^ a^), the scalar part of the pseudopotential de- 
termines the "band renormalization" of l/Ti, that is, the 
enhancement due to spin hot spots . Here we can only 
offer a guess. Considering the spin-orbit part "fixed," 
our semi-empirical scalar pseudopotential, which repro- 
duces the experimental band gaps at symmetry points 
within 5%, does not introduce more than another 10% 
error |Q , making A5 determined with 20% accuracy. As 
for the spin-orbit interaction, future experiments done 
in the regime where Ti ~ 1/T (that is, above 200 K), 
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FIG. 2. Spin-flip Eliashberg function OgF for aluminum. 
The dotted line shows the phonon density of states F and the 
dashed line is the ordinary (spin-conserving) Eliashberg func- 
tion a^F. The curves are calculated for the model described 
in the text. 

will have the opportunity to set definite constraints on a 
through a direct comparison with our theory. 

In summary, we have provided the first fully quantita- 
tive theory for the temperature dependent spin relaxation 
rate in aluminum taking into account spin-orbit coupling 
and electron-phonon interaction within the Elliott- Yafet 
formalism using realistic pseudopotentials. Our theoret- 
ical results are in excellent agreement with the measured 
Ti{T) in aluminum and for T > 100 K, where experimen- 
tal results are currently non-existent, our theory provides 
specific predictions for comparison with future experi- 
ments. 
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